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INTRODUCTION 


This  is  a  final  report  on  collaborative  study  between  AeroChem  Research  Laboratories 
(principal  investigator  —  Dr.  H.  F.  Calcote)  and  the  Pennsylvania  State  University  (principal 
investigator  —  Dr.  M.  Frenklach).  The  ultimate  objective  of  this  program  was  to  develop  a 
quantitative  physical/chemical  model  of  the  mechanism  of  soot  formation  that  is  consistent  with 
available  experimental  data.  The  specific  objectives  of  the  proposed  3-year  study  were:  1)  to 
delineate  the  relative  importance  of  radical  and  ionic  mechanisms  of  soot  nucleation  in  flames, 
and  2)  to  determine  the  optimum  chemical  and  physical  model  of  the  total  soot  formation 
process.  Only  the  work  performed  at  Penn  State  is  reported  here,  that  of  AeroChem  is  reported 
in  their  parallel  final  technical  report 

WORK  PROPOSED 

The  general  objectives  proposed  for  the  Penn  State  efforts  were: 

1.  To  undertake  flame  simulation  and  reaction-path  analysis  using  the  ion-reaction  mechanism 
provided  to  us  by  AeroChem. 

2.  To  analyze  the  competition  for  the  formation  of  soot  and  soot  precursors  between  the  ionic 
and  neutral  reaction  mechanisms. 

3.  To  update  the  Penn  State  neutral-reaction  mechanism  for  the  formation  of  polycyclic 
aromatic  hydrocarbons  (PAH). 

4.  To  develop  a  numerical  algorithm  and  a  computer  code  for  modeling  the  formation  and 
growth  of  large-sized  species  and  soot  particles. 

5.  To  perform  computer  simulation  of  soot  particle  formation  in  flames  with  the  code 
developed. 


RESULTS 

During  the  three-year  period  of  the  project,  we  completed  most  of  the  planned  work: 
several  ion  reaction  mechanisms  provided  to  us  by  AeroChem  were  “cleaned-up”  and  run  with  a 
flame  code,  the  computational  results  were  analyzed,  a  computer  code  for  modeling  large 
species  and  soot  particle  formation  and  growth  was  developed,  our  neutral  reaction  mechanism 
was  updated  and  tested,  and  a  detailed  computer  simulation  of  soot  particle  nucleation  and 
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growth  in  laminar  premixed  flames  was  successfully  performed  using  our  updated  neutral 
reaction  mechanism.  The  latter  part  is  a  particularly  significant  accomplishment,  as  it  is  for  the 
first  time  soot  particle  inception  could  be  modeled  from  first  principles,  starting  with  fuel 
decomposition  and  going  all  the  way  to  predicting,  in  quantitative  agreement  with  experiment, 
soot  particle  properties.  The  following  is  a  detailed  account  for  the  work  performed  and  the 
results  attained. 

Ion  Reaction  Mechanism 

We  began  our  work  by  testing  the  general  performance  of  an  ion-molecule  reaction 
mechanism  under  shock-tube  conditions.  The  premise  of  such  testing  is  founded  on  our 
previous  result  that  the  major  neutral-reaction  pathways  for  PAH  formation  under  shock-tube 
oxidation1  and  laminar  premixed  flame2  conditions  are  very  similar  to  each  other.  The  initial 
ionic  mechanism  provided  to  us  by  AeroChem  has  been  reported  and  discussed  in  the 
AeroChem  reports.3*4  The  detailed  account  of  the  shock-tube  calculations  is  given  in  our  first 
annual  report.5 

The  basic  results  of  the  shock-tube  simulations  indicated  that  the  ionic  mechanism  produces 
polycyclic  aromatic  hydrocarbons  at  significantly  lower  rate  than  does  the  mechanism  composed 
of  only  reactions  of  neutral  species.  An  extensive  sensitivity  and  rate  analysis  found  no 
conceivable  adjustments  in  the  parameter  values  that  would  affect  this  conclusion.  There  are 
two  basic  factors  identified  by  the  sensitivity  and  rate  analysis  to  be  responsible  for  the  low  rate 
of  PAH  production  via  the  ionic  mechanism  tested.  First  is  the  relatively  low  rate  of  ion 
production.  The  second  factor  limiting  the  rate  of  PAH  production  via  the  ionic  mechanism  is 
the  reversibility  of  some  of  the  principal  reaction  steps,  and  a  significant  number  of  steps  in  the 
principal  reaction  sequence  was  found  to  be  reversible.  The  computer  simulations  indicated  that 
the  forward  and  reverse  reactions  of  these  steps  are  tightly  balanced,  Lc.,  these  reactions  are  in 
partial  equilibrium,  similarly  to  the  case  for  the  neutral-reaction  mechanism.2*6  This  reduces 
dramatically  the  net  mass  flux  and  thus  limits  the  rate  of  PAH  production.  Among  other  things, 
it  implies  that  an  oversimplified  analysis  presented  recently,7  ignoring  the  kinetic-thermodynamic 
coupling  controlling  the  rate  of  the  principal  reaction  pathway,  cannot  possibly  be  valid. 

We  next  examined  the  competition  between  the  ionic  and  neutral  mechanisms  for  PAH 
production  at  flame  conditions.  For  this  purpose,  a  “well-studied"  laminar  premixed  flame,  for 
which  the  experimental  data  on  ionic  species  concentrations  are  available,  had  to  be  chosen  for 
the  computer  simulations.  As  there  is  no  single  flame  reported  for  which  there  are 
measurements  available  of  all  the  pertinent  chemical  species,  two  close  flames  were  chosen  to 
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represent  a  “well-studied”  flame.  As  the  primary  focus  of  this  study  is  on  the  ionic  species,  the 
main  flame  —  referred  to  as  Flame  II  —  is  the  flame  for  which  ion  concentration  profiles  were 
measured  by  Calcote  and  co-workers:8  52.9  %  C2H2  -  44.1  %  02-Ar  (0  =  3.0),  pressure  20 
torr,  cold  gas  velocity  50  cm/s.  Since,  however,  Calcote  and  co-workers  did  not  measure  the 
concentrations  of  neutral  species  in  this  flame,  a  second  flame,  that  of  Delfau  and  Vovelle9  — 
referred  to  as  Flame  12  —  was  chosen  to  test  the  predictions  of  the  present  simulations  for  the 
major  neutral  species.  The  stoichiometry,  pressure  and  cold  gas  velocity  of  this  flame  are 
exactly  the  same  as  those  of  Flame  II.  The  mixture  used  by  Delfau  and  Vovelle  did  not  contain 
any  argon,  whereas  that  of  Calcote  and  Keil  had  3  %  argon.  Also,  the  reported  temperature 
profiles  of  Flames  II  and  12  are  slightly  different  from  each  other.  Temperature  profile  is  a 
critical  input  parameter  for  the  flame  simulations  and  even  small  differences  in  the  given  profile 
data  cause  significant  changes  in  computed  profiles  of  minor  (e.g.,  PAHs)  species. 

The  flame  computations  were  performed  using  the  Sandia  burner  code.10  Several  reaction 
mechanisms  provided  to  us  by  AeroChem  have  been  tested  and  analyzed  during  the  course  of 
the  study.  Some  of  the  results  of  these  flame  simulations,  thsoe  obtained  with  the  most  recent 
version  of  the  AeroChem  mechanism,  are  presented  in  Figures  1  through  4.  Examination  of 
Figure  1  indicates  that  the  concentrations  of  major  species  are  predicted  well  —  as  good  as  the 
present  state  of  the  art  in  the  knowledge  of  chemical  reactions,  their  rates  and  thermochemistr  , 
the  accuracy  of  numerical  techniques,  and  quality  of  experimental  measurements  allow.  The 
comparison  between  the  computed  and  experimental  concentrations  of  ions  is  given  in  Figures  2 
and  3.  The  agreement  is  worse  for  these  species  compared  to  the  agreement  seen  for  the  major 
species  in  Figure  1.  In  the  computer  simulation  as  opposed  to  the  experiment,  the  positions  of 
the  concentration  peaks  are  shifted  towards  the  burner  and  the  decay  in  the  concentration 
profiles  is  much  steeper.  Similar  behavior  was  seen  in  computer  simulations  of  an  acetylene 
flame  with  the  neutral-reaction  mechanism.2  The  reasons  for  this  are  discussed  in  the  work 
cited.  It  is  important  to  notice,  however,  that  the  peak  values  of  the  ion  concentrations  are 
predicted  reasonably  well,  taking  into  account  the  uncertainty  in  the  input  reaction  data  and  in 
the  experimental  measurements  —  most  of  them  are  computed  within  an  order  of  magnitude  of 
the  experimental  numbers.  In  this  sense,  the  results  obtained  with  different  versions  of  the 
AeroChem  ionic  mechanism,  as  far  as  the  ions-vs-neutrals  competition  is  concerned,  are  very 
similar  to  each  other  (compare,  for  example,  the  data  in  present  Figures  1  through  4  with  the 
data  reported  in  Figures  5  through  8  in  our  first  annual  report5). 

Although  the  analysis  of  the  results  obtained  in  the  flame  simulations  is  not  entirely 
completed,  the  following  conclusion,  similar  to  to  those  reached  based  on  the  shock-tube 
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figure  2.  Concentration  profiles  of  small  ions  in  Flame  11:  a)  experimental  measurements;8  (b)  present  computations  using  the 
temperature  profile  of  Ref.  8. 
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Profiles  of  selected  aromatic  hydrocarbons  computed  with  radical  (solid  lines)  and  ionic  (dotted  lines) 
mechanism. 


simulation  results,  has  emerged:  the  ionic  mechanism  produces  polycyclic  aromatic 
hydrocarbons  at  significantly  lower  rates  than  does  the  neutral-reaction  mechanism.  This  result 
is  illustrated  in  Fig.4,  where  the  concentrations  of  selected  aromatics  computed  with  the  ionic 
mechanism  are  compared  to  those  computed  with  the  neutral-reaction  mechanism  at  the  same 
flame  conditions.  As  can  be  seen  from  this  comparison,  the  production  of  aromatics  via  the 
neutrals  is  much  faster  than  that  via  ions  (the  difference  between  the  two  cases  increases  with 
the  switching  off  oxidation  reactions  in  computer  simulations  with  the  neutral-reaction 
mechanism;  also  note  that  the  computed  concentrations  of  ions  are  significantly  larger  than  the 
experimental  counterparts).  A  sensitivity  and  rate  analysis  performed  so  far  found  no 
conceivable  adjustments  in  the  parameter  values  that  would  change  this  principal  conclusion. 
The  factors  responsible  for  this  are  the  same  as  those  identified  in  the  analysis  of  the  shock-tube 
results:  a  relatively  low  rate  of  chemi-ionization  as  compared  to  that  of  radical  production  at  the 
same  conditions,  and  reversibility  of  the  reaction  steps  responsible  for  the  formation  of  high- 
molecular  weight  species. 

Another  shortcoming  of  the  present  ionic  mechanism  is  the  lack  of  knowledge  on  the 
microscopic  mechanism  of  many  ion-molecule  reaction  in  the  AeroChem  scheme.  In  fact,  it  is 
unlikely  that  many  individual  reaction  steps  will  proceed  as  written  and  with  the  rate  coefficient 
values  assigned  to  them.  Most  of  these  ion  reactions  are  written  on  the  premises11  that  ion- 
molecule  reactions  have  no  potential-energy  barriers  and  that  ionic  species  isomerize  fast 
However,  some  available  (although  limited)  infotmation  indicates  to  the  opposite.  For  instance, 
we  performed  quantum-chemical  calculations  on  a  key  ion-molecule  step.  It  appears  that  even  a 
more  conceivable  reaction  step  of  the  AeroChem  mechanism  has  a  measurable  potential  energy 
barrier.  This  means,  that  the  rate  coefficients  and,  more  importantly,  reaction  channels  of  the 
present  ionic  mechanism  are  most  likely  to  greatly  overestimate  the  growth  rate  of  polycyclic 
aromatic  species  in  flame  environments. 

Neutral  Reaction  Mechanism 

Accurate  prediction  of  polycyclic  aromatic  hydrocarbons  in  flames  is  critical  for  modeling  of 
soot  formation  in  hydrocarbon  combustion.  In  our  initial  efforts,  the  objective  was  a  proof-of- 
a-concept  computer  simulation  of  PAH  production  in  high-temperature  hydrocarbon  pyrolysis,6 
and  later,  in  a  low-pressure  laminar  premixed  acetylene  flame.2  Now,  when  the  acetylene- 
addition  mechanism  we  have  advocated  for  the  last  seven  years  is  well  tested  under  various 
experimental  conditions  and  by  different  research  groups,  the  objective  is  to  improve  the 
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mechanism  to  a  point  that  it  is  capable  of  quantitative  predictions.  Initial  results  of  such  efforts 
are  reported  below. 

We  chose  to  model  the  PAH  profiles  in  three  well-studied  laminar  premixed  flames,  whose 
conditions  are  summarized  in  the  following  table: 


Table  1.  Flame  conditions  for  modeling  with  neutral-species  mechanism 


No. 

Composition  in  Argon 

P  (torr) 

Cold  gas  velocity  (cm/s) 

Reference 

N1 

46.5%  C2H2  -  48.5%  O2 

20 

50 

MIT  group12 

N2 

23.6%  C2H2  -  21.4%  O2 

90 

20 

Bockhom  et  al. 13 

N3 

16.5%  C2H4  -  17.9%  O7 

760 

7.8 

Harris  etalM 

These  flames  were  computed  using  the  Sandia  burner  code.10  The  Lennard- Jones  parameters  of 
PAH  molecules  were  calculated  using  empirical  correlations  and  available  physical  properties.15 
The  transport  coefficients  obtained  in  this  manner  agreed  well  with  experimental  data  available 
for  several  compounds.  The  thermodynamic  data  were  taken  from  Stein  and  Fahr,16  except  for 
the  Q-H  group  additivity  which  was  based  on  the  recent  results  of  Green  et  a/.17 

The  reaction  mechanism  is  composed  of  two  parts:  the  small  molecule  reactions,  responsible 
for  the  main  flame  structures;  and  reactions  describing  the  formation  and  growth  of  PAHs.  The 
reaction  mechanism  for  the  former  part  was  composed  of  data  taken  primarily  from  several 
sources:  (a)  our  recently  optimized  methane  oxidation  mechanism,18  (b)  revised  mechanism  of 
Frenklach  and  Wamatz,2  and  (c)  revised  acetylene  oxidation  mechanism  of  Hwang  era/.19  The 
main  revisions  included:  the  rate  coefficients  of  acetylene-hydroxyl,  and  acetylene-oxygen  atom 
reactions 

C2H2  +  OH  products 

C2H2  +  O  -»  products 

were  taken  from  Miller  and  Melius,20  and  Mahmud  and  Fontijn,21  respectively.  The  rate 
coefficient  for  reaction 

C2H2  +  C2H3  C4H4  +  H 
was  taken  from  Fahr  and  Stein.22 

The  formation  of  the  first  aromatic  ring  followed  the  basic  reaction  scheme  of  Frenklach  and 
Wamatz,2 

n-C4H3  +  C2H2  ->  phenyl 
n-GtHs  +  C2H2  -»  benzene  +  H 
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however  with  the  rate  coefficients  taken  from  the  work  of  Westmoreland  et  ai.,23  who  used  a 
QRRK  methodology  of  estimating  chemical  activation  processes  involved  in  these  reactions. 
The  assignment  of  the  rate  coefficients  for  the  main  reaction  classes  responsible  for  PAH 
growth  are  summarized  in  Table  2. 


Table  2.  Revision  of  rate  coefficients  for  main  reaction  classes  of  PAH  growth 


Reaction 

Prototype  reaction 

Reference 

aryl-H  +  H  -4  aryl  +  H2 

benzene  +  H  -4  phenyl  +  H2 

Kiefer  ef  a/.24 

aryl  +  C2H2  -4  aryl-C2H  +  H 

phenyl  +  C2H2  -4  phenyl-acetylene  +  H 

Fahr  and  Stein22 

aryl  +  O2  products 

phenyl  +  O2  -4  products 

Lin  and  Lin25 

The  computational  results  for  the  three  flames  are  presented  in  Figures  5  to  10.  These 
results  were  obtained  with  the  rate  coefficients  for  several  small-molecule  reactions  adjusted  by  a 
factor  of  2  to  3,  still  within  their  experimental  uncertainties.  As  can  be  seen  in  these  figures, 
the  agreement  between  the  modeling  predictions  and  experimental  data  is  good  for  major 
species.  The  agreement  for  PAHs  is  not  exactly  quantitative,  but  certainly  much  closer  —  in 
both  shape  and  absolute  values  —  than  reported  by  us  previously.2  The  important  result  here  is 
that  the  near-quantitative  accuracy  was  obtained  for  PAHs  with  a  single  reaction  mechanism  for 
all  the  three  flames  simultaneously.  The  major  reaction  pathways,  their  roles  in  PAH  growth, 
and  the  sensitivity  information  are  essentially  the  same  as  reported  and  discussed  in  the  previous 
work.2 
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Figure  5.  Computed  vs.  experimental  major  Figure  6.  Computed  vs.  experimental  major 
species  concentrations  of  Flame  Nl.  aromatic  species  concentrations  of  Flame  Nl. 
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Mole  fraction 


0.0  0.5  1.0  1.5  2.0  2.5 

x(cm) 


Figure  7.  Measured  mole  fractions  of  diacetylene 
and  selected  PAHs  of  Flame  N2. 
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Figure  8.  Computed  mole  fractions  of 
diacetylene  and  PAHs  of  Flame  N2. 


0.0  0.1  02  0.3  0.4  0.0  0.1  0.2  0.3  0.4 

x  (cm)  x  (cm) 

Figure  9.  Computed  vs.  experimental  major  Figure  10.  Computed  vs.  experimental 

species  concentrations  of  Flame  N3.  minor  species  concentration  of  Flame  N3. 
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Flame  Simulation  of  Soot  Particle  Nucleatlon  and  Growth 


Using  the  neutral  reaction  mechanism  described  in  the  previous  section,  we  performed 
detailed  modeling  of  soot  particle  nucleation  and  growth  in  several  laminar  premixed  flames. 
The  computational  model  consists  of  three  logical  parts:  (I)  initial  PAH  formation ,  which 
includes  a  detailed  chemical  kinetic  description  of  acetylene  pyrolysis  and  oxidation,  formation 
of  the  first  aromatic  ring,  and  its  subsequent  growth  to  a  prescribed  size;  (II)  planar  PAH 
growth,  comprised  of  replicating-type  growth  of  PAHs  beyond  the  prescribed  size;  and  (IH) 
spherical  particle  formation  and  growth,  consisting  of  coagulation  of  PAHs  formed  in  part  (II) 
followed  by  the  growth  of  the  resulting  particles  by  coagulation  and  surface  reactions. 

The  initial  part  of  the  particle  dynamics  model,  nucleation,  describes  the  planar  growth  of 
PAHs  via  the  H-abstraction/C2H2-addition  reaction  sequence2*6  using  the  technique  of  chemical 
lumping.26*27  This  method  provides  a  mathematically  rigorous  description  of  the  growth 
process  to  an  infinite  size  PAH,  which  can  be  schematically  represented  as 

A /  -*  A i+i  -»  A i+2  ->  ->  Aoo  ,  (1) 

where  A/  represents  an  aromatic  species  containing  1  fused  rings.  The  PAH  species  formed  in 
(1)  are  then  allowed  to  coagulate,  that  is,  all  the  A/s  (i  =1,  ffl„..,<»)  collide  with  each  other 
forming  dimers;  the  dimers,  in  turn,  collide  with  Aj  forming  trimers  or  with  other  dimers 
forming  tetramers;  and  so  on.  The  coalescence  reactions  were  treated  as  irreversible  having 
sticking  coefficients  of  unity.  As  the  focus  of  this  work  was  on  very  young,  small  particles,  it 
*vas  assumed  that  the  coagulation  dynamics  is  in  the  free-molecular  regime.  A  size-independent 
enhancement  factor  of  2.2,  based  on  the  results  of  Hams  and  Kennedy,28  was  used  in 
calculations  of  collision  frequencies. 

Beginning  with  the  dimers,  the  forming  clusters  were  assumed  to  be  “solid  phase”  and 
allowed  to  add  and  lose  mass  by  surface  reactions 


Csoot-H 

+  H  < 

Csoot*  +  H2 

(2) 

Qoot*  + 

H  -> 

Csoot-H 

(3) 

Csoot*  + 

c2h2 

->  Csoot-H  +  H 

(4) 

Csoot*  + 

02  -» 

products 

(5) 

Csoot-H 

+  OH 

-»  products  , 

(6) 

where  Cjoqj-H  represents  an  arm-chair  site  on  the  soot  particle  surface  and  C^t*  the 
corresponding  radical.  This  mechanism  is  adopted  based  on  the  postulate29-31  that  the 
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H-abstraction/C2H2-addition  (HACA)  reaction  sequence6  is  responsible  for  high-temperature 
growth  of  all  forms  of  carbonaceous  materials.  Wiese hnowsky  et  al?2  also  suggested  that  the 
HACA  reactions2*6  “offer  a  key  to  the  understanding  of  a  number  of  phenomena”  observed  in 
their  flame  study. 

Estimation  of  the  rate  coefficients  of  the  heterogeneous  reactions  was  based  on  analogous 
gas-phase  reactions  of  one-ring  aromatics,  benzene  and  phenyl.  The  reference  gaseous  reactions 
for  heterogeneous  reactions  (2)-(5)  were  the  corresponding  gas-phase  reactions  in  Table  (2). 
Oxidation  by  OH,  reaction  (6),  was  assumed  to  proceed  with  the  collision  efficiency  of  0.13, 
based  on  the  results  of  Neoh  et  al.33  The  particle  dynamics  —  the  evolution  of  soot  particles 
undergoing  simultaneous  nucleation,  coagulation  and  surface  reactions  described  above  —  was 
modeled  by  a  method  of  moments  (Method  II  of  Ref.  34).  This  method  does  not  require  the 
assumption  of  a  particle  size  distribution  function  (PSDF).  The  closure  of  the  differential 
equations  for  the  PSDF  moments  is  accomplished  by  interpolation  between  the  moments.26*34 
The  numerical  integration  of  the  moment  equations  is  extremely  economical;  for  instance,  one 
computer  run  took  from  20  to  50  s  on  an  IBM  3090/600S  main-frame  computer. 

Figures  11  and  12  present  the  results  computed  using  the  model  discussed  above  for  two  of 
tite  chosen  flames:  Flame  B  —  flame  No.  1  of  Bockhom  and  co-workers32 — 25.4  %  C2H2  - 
19.6  %  O2  —  Argon,  pressure  90  torr,  and  Flame  H  —  an  atmospheric  16.5%  C2H4  -  17.9  % 
O2  -  Argon,  of  Harris  and  co-workers.35  Both  flames  woe  computed  with  the  same  kinetic 
model,  assuming  A4  (pyrene  or  acepyrene)  to  be  the  A/  species  which  initiates  the  planar 
growth  of  PAHs  and  the  formation  of  particles,  parts  (II)  and  (ID)  of  our  model,  respectively. 
Changing  A/  affects  strongly  the  rate  of  particle  nucleation. 

It  can  be  seen  in  Figs.  1 1  and  12  that  the  model  predictions  are  in  relatively  close  agreement 
with  experiment  for  the  initial,  particle  inception  part  of  the  flames;  however,  the  computed  and 
experimental  results  begin  to  deviate  from  one  another  with  increasing  flame  height  One  factor 
identified  as  a  cause  of  this  disagreement  is  the  shape  of  the  PAH  profiles.  This  is 
demonstrated  in  Fig.  11,  where  it  is  shown  that  using  a  narrower  (experimental-like13)  profile 
for  A4  improves  the  agreement  between  the  model  and  experiment  for  Flame  B.  It  was  also 
computed  that  PAH  condensation  on  the  particle  surface  is  sufficiently  fast  to  significantly  affect 
the  PAH  profiles  in  the  post-flame  zone. 
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Fig.  11.  Comparison  of  model  predictions  and  experimental  data  far  Flame  B:  (a)  particle  number  density; 
(b)  particle  surface  area;  (c)  average  particle  diameter  [circles  -  experimental  results,32  solid  lines  -  model  prediction, 
dotted  lines  -  model  prediction  using  an  experimental-like  profile  for  A4];  (d)  specific  rate  of  soot  panicle  oxidation 
by  02  [solid  line  -  model  prediction,  dashed  line  -  from  the  expression  of  Nagle  and  S  trickland-Cons table.36 
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Fig.  12.  Comparison  of  model  predictions  and  experimental  data  for  Flame  H:  (a)-(c)  N  -  particle  number  density, 
r  -  particle  radius,  <rp>  -  average  d  -  average  particle  diameter,  circles  -  experimental  results,' 55  solid  lines  - 
model  prediction;  (d)  specific  surface  growth  rate  coefficient. 
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As  expected  from  the  dynamics  of  Smoluchowsld  coagulation,  the  computed  rate  of 
nucleation  is  balanced  by  the  rate  of  coagulation  throughout  the  particle  inception  zone.  Particle 
inception  is  primarily  determined  by  PAH  coagulation,  initiated  and  controlled  by  PAH 
coalescence  into  dimers.  Our  model  describes  the  surface  processes  in  terms  of  elementary 
reactions  of  active  sites.  The  density  of  active  sites  is  determined  by  the  gas-phase  environment 
—  by  the  competition  among  reactions  (2)-  (6).  Hence,  the  dissimilarity  in  surface  growth 
behavior  can  be  explained  by  the  difference  in  the  dynamics  of  active  sites.  Our  model  resolves 
the  controversy  brought  up  recently  by  Bockhom  and  co- workers,32  who  observed  that  the 
surface  growth  rate  is  independent  of  the  surface  area  of  soot  particles,  contrary  to  Harris  and 
Weiner,37  who  postulated  that  the  rate  of  surface  growth  is  first  order  in  particle  surface  area. 
Also,  the  computed  net  surface  growth  rate  is  computed  to  be  in  close  agreement  with  that 
determined  by  Harris  and  Weiner32  (Fig.  12d)  and  the  predicted  rate  of  soot  oxidation  by  O2 
agrees  well  with  the  expression  of  Nagle  and  S trickland-Cons  table36  (Fig.  lid). 

Our  model  predicts  what  is  essentially  the  classical  picture  of  soot  particle  inception:38 
formation  of  PAHs,  their  coagulation  into  clusters  and  cluster  coagulation,  and  surface  growth 
due  to  chemical  and  physical  processes.  The  model  also  predicts  the  classical  structure  of  soot 
particles:39  a  less  dense  particle  care,  composed  of  randomly  oriented  PAH  oligomers,  and  a 
more  dense  concentrically-arranged  particle  shell. 

PRINCIPAL  CONCLUSIONS  AND  RESULTS 

The  principal  conclusion  of  this  study  is  that  the  formation  of  polycyclic  aromatic  precursors 
to  soot  in  hydrocarbon  flames  is  dominated  by  the  reactions  of  neutral  species  and  not  ionic 
ones.  Whether  ions  play  role  in  other  soot-forming  processes,  like  PAH  coagulation  is  yet  to 
be  investigated. 

A  computer  code  for  modeling  large  species  and  soot  particle  formation  and  growth  was 
developed. 

A  detailed  computer  simulation  of  soot  particle  nucleation  and  growth  in  laminar  premixed 
flames  was  successfully  performed  using  our  updated  neutral  reaction  mechanism  —  for  the 
first  time  soot  particle  inception  could  be  modeled  from  first  principles,  starting  with  fuel 
decomposition  and  going  all  the  way  to  predicting,  in  quantitative  agreement  with  experiment, 
soot  particle  properties. 
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FUTURE  WORK  DIRECTION 


Following  the  successful  accomplishments  in  soot  formation  model  development,  culminated 
in  quantitative  prediction  of  soot  particle  properties  in  several  laminar  premixed  flames,  the 
direction  for  the  future  research,  under  the  new  AFOSR-  sponsored  project  entitled 
“Development  of  Predictive  Reaction  Models  of  Soot  Formation,”  is  to  extend  the  modeling 
efforts  to  computer  simulation  of  soot  particle  phenomena  in  other  flames  of  hydrocarbon  fuels 
along  with  further  refinement  of  the  underlying  reaction  mechanism  of  soot  formation. 
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